Particle spectra in Pb-Pb collisions at ^/snn = 2.76 TeV 
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Particle production in Pb-Pb collisions at ^snn ~ 2.76TeV is studied in the (3 + l)-dimensional 
viscous hydrodynamic model. The shapes of the calculated transverse momentum spectra of tt"^, 
K'^ , protons, H~, and Q~ are in satisfactory agreement with preliminary data of the ALICE Col- 
laboration, while the particle ratio p/vr^ is shghtly overpredicted, and the strange barion yields are 
underpredicted . 
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I. INTRODUCTION 

Particle production in Pb-Pb collisions at the high- 
est available energy ^/snn = 2.76 TeV has been stud- 
ied experimentally at the CERN Large Hadron Collider 
(LHC) [1-3]. The observation of the elliptic and triangu- 
lar flows indicates that a collectively expanding fireball 
of dense matter is formed, confirming results obtained in 
collisions at lower energies [4] . The hydrodynamic model 
of the dynamics provides a quantitative explanation for 
observables defined for particles emitted with soft mo- 
menta [5]. In particular, hydrodynamic models are ap- 
plied to describe the anisotropic flow of charged particles 
produced in Pb-Pb collisions at the LHC [6-14]. 

Statistical models of the particle production in heavy- 
ion collisions predict the production rates of specific 
hadrons assuming a chemically equilibrated system [15]. 
Recent results for Pb-Pb collisions at the LHC seem to 
be incompatible with this simple mechanism. The num- 
ber of protons emitted is lower than expected from the 
rates of the emission of other particles, which raises a 
doubt in the assumption that the production of different 
particle species happens at a common chemical freeze- 
out temperature [16]. Particle abundances may undergo 
significant modifications in the nonequilibrium dynam- 
ics after hadronization, e.g. annihilation processes may 
reduce baryon multiplicities [17]. 

Transverse momentum spectra of identified particles 
represent a more basic observable, as they involve the 
overall multiplicity as well as the momentum distribu- 
tion for each particle species. The hydrodynamic evo- 
lution of a fireball of thermally equilibrated fluid until 
the freeze-out temperature Tf, leads to a common chem- 
ical and kinetic freeze-out temperature. A similar idea 
constitutes the basic assumption of the single freeze-out 
model of particle emission [18]. Additional freedom is al- 
lowed in hydrodynamic models assuming that below the 
chemical freeze-out temperature the matter is in kinetic 
equilibrium but particle abundances remain frozen [19]. 
In such models particle ratios correspond to a fixed chem- 
ical freeze-out temperature Tchem, while the transverse 
momentum spectra are determined by the convolution 



of the collective velocity with the thermal emission at 
the kinetic freeze-out temperature Tkm- The genuinely 
nonequilibrium phase at the end of the evolution may 
be addressed using a hybrid model with hydrodynamics 
for the dense phase and a hadronic cascade afterburner 
for the latter evolution [20-24]. The evolution in the 
hadronic cascade changes the chemical composition in 
the system and the momentum spectra of particles. 



Deviations from local equilibrium in the hydrodynamic 
model of the dynamics are introduced as viscosity correc- 
tions to the energy- momentum tensor [25-35] . In partic- 
ular, shear viscosity is important in quantitative predic- 
tions for the elliptic and triangular collective flows. The 
rapid expansion of the fireball introduces sizable correc- 
tions from bulk viscosity, if the equilibration processes 
are not fast enough to restore local equilibrium. Such 
deviation are twofold. First, the chemical composition 
remains effectively frozen at some stage, while the local 
energy density drops. Second, the local momentum dis- 
tributions of particles in the fluid become softer. Bulk 
viscosity leads to both effects, depending on the local 
expansion rate [34, 36-38] . 



We present calculations for the transverse momentum 
spectra of identified particles in a (3 -I- I)-dimensional 
[(3 -I- 1)-D] hydrodynamic model with bulk and shear 
viscosities for Pb-Pb collisions at LHC energies. The 
presence of bulk viscosity in the hadronic phase yields 
nonequilibrium effects for the chemical composition and 
for the transverse momentum spectra. We find a satis- 
factory agreement with preliminary data of the ALICE 
Collaboration [3, 16] for pion, kaon, proton spectra and 
abundances. At the same time we constrain the freeze- 
out temperature using the interferometry data [2]. Re- 
producing the charged particle density in pseudorapidity 
[3] and the transverse momentum spectra allows for a di- 
rect prediction on the rapidity distributions of identified 
hadrons. 



II. INITIAL CONDITIONS AND 
HYDRODYNAMIC EVOLUTION 

The expansion of the firebah is described using sec- 
ond order viscous hydrodynamics. Hydrodynamic equa- 
tions are solved in (3 + 1)-D together with the Israel- 
Stewart equations for the stress corrections tt^*^ and 11 to 
the energy-momentum tensor (for details see [37]). We 
use a constant ratio of shear viscosity to entropy den- 
sity r]/s — 0.08. The bulk viscosity is nonzero only in 
the hadronic phase, we use (/s — 0.04 and C/s — 0.08. 
The equation of state is a combination of the lattice QCD 
[39] and hadron gas equations of state, obtained in a ther- 
modynamically consistent way [40]. The initial time for 
the hydrodynamic expansion is 0.6 fm/c. The relaxation 
times in the Israel-Stewart equations are r,r = Tn = -j^. 
For the chosen values of the relaxation times the stress 
corrections are close to the ones from the Navier-Stokes 
expression, at latter stages of the expansion. For a strong 
temperature dependence of the viscosity coefhcients, the 
value of the relaxations time can influence the evolution 
[41]. 

The initial entropy density s{rn^,x,y) for the (3 -I- 1)- 
D hydrodynamic evolution in the space-time rapidity rjii 
and the transverse coordinates a;, y is 
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In the transverse plane the density is defined as a 
combination of the participant nucleon ppart = -^-i- + 

densities. The factor 

— (jv +jv *) 1 ^^ ^'i' ^ implements in the ini- 
tial density the assumption that forward (backward) go- 
ing nucleons N^ (iV_ ) emit particles preferentially in the 
forward (backward) rapidity hemisphere [42]. The pa- 
rameters of the longitudinal profile 
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are adjusted to reproduce the charge particle density in 
pseudorapidity, 770 = 2.3 and ct,, = 1.4, yi, is the beam 
rapidity. The parameters of the Glauber model used to 
calculate the entropy profiles for Pb-Pb collisions {A = 
208) are Ra = 6.48fm and a — 0.535fm; the nucleon- 
nucleon cross section is cr = 62mb. 

At the freeze-out temperature particles are emit- 
ted from the freeze-out hypersurface according to the 
Cooper-Frye formula with viscosity corrections. The 
nonequilibrium modifications of the equilibrium momen- 
tum distribution /o from shear viscosity are quadratic in 
momentum [26] 
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FIG. 1. (Color online) Charge particle distributions in pseu- 
drapidity in Pb-Pb collision at \/sjv]v ~ 2.76 TeV, for cen- 
tralities (from top to bottom) - 5%, 5 - 10%, 10 - 20%, . . . , 
50 — 60%. Results from viscous hydrodynamics (lines) are 
compared to preliminary data of the ALICE Collaboration 
(symbols) [3]. 



The correction from bulk viscosity are taken from the 
relaxation time formula [34, 43] 
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where c^ is the sound velocity and Chuik is a normaliza- 
tion constant. At the freeze-out hypersurface, the flow 
velocity and the stress corrections from viscosity are ex- 
ported to a Monte-Carlo statistical emission code [44] 
that is used to generate particle spectra. Bulk viscosity 
corrections are calculated with respect to the equilibrium 
distribution having the same energy density. Imposing 
constraints corresponding to conserved charges leads to 
deviations from the chemical equilibrium with respect 
to this equilibrium reference state. Effectively, it means 
that due to incomplete equilibration the chemical equilib- 
rium temperature is shifted up and the effective kinetic 
temperature is shifted down (due to a redshifting of par- 
ticle thermal momenta). An expanded discussion of this 
approach can be found in [37, 38]. 



III. RESULTS 

The parameters of the initial entropy density for the 
hydrodynamic evolution are constrained by the charge 
particle density as function of pseudorapidity measured 
for different centrality classes (Fig. 1). We find that the 
optimal mixing parameter for the admixture of binary 
collisions a — 0.15; the same value as used in (2 + 1)-D 
viscous model simulations for heavy-ion collisions at the 
LHC [10]. The half-width of the plateau in the distribu- 
tion rjQ = 2.3 is larger as compared to tjq = 1.5 that has 
been used for Au-Au collisions at ^/sWn = 200 GeV [37]. 
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FIG. 2. (Color online) Transverse momentum spectra of 
■k'^ , K^ , p (from top to bottom) in Pb-Pb collisions with 
centralities — 5%, obtained in the viscous hydrodynamic 
model with (/s = 0.04 and T/ = 150 MeV (dashed lines), 
C/s = 0.04 and Tf = 140 MeV (dotted lines), C/s = 0.08 and 
Tf — 140 MeV (solid lines). Symbols represent preliminary 
data of the ALICE Collaboration [3] . The dashed-doted line 
represents the result of the viscous hydrodynamic calculation 
with C/s = 0.08 and Tf = 140 MeV but without bulk viscosity 
corrections at freeze-out. 



Transverse niomentum spectra for p"^, K'^ and pro- 
tons in central collisions are presented in Fig. 2. The 
spectra get harder when the freeze-out temperature is 
lowered from Tf = 150 MeV (dashed lines) to 140 MeV 
(dotted lines), and using (^/s — 0.04. Especially, the pion 
spectra get too flat. The size of the fireball at freeze-out 
can be estimated from the interferometry radii [45] . The 
radius Rside is best described using a freeze-out temper- 
ature of 140 MeV (Fig. 3). Also the value of the ratio 
Rout/ Rside is closcr to the experimental data if the evo- 
lution is longer. The proton yield is too large for the 
freeze-out at 150 MeV, while the average transverse mo- 
mentum of protons is too low. The proton spectra and 
yields are better reproduced for Tf — 140 MeV. To de- 
scribe at the same time the pion and proton spectra a 
freeze-out temperature 140 MeV and a bulk viscosity co- 
efficient C/s = 0.08 are used (solid lines in Figs. 2 and 3). 
Increasing the bulk viscosity makes the local momentum 
distributions of light particles (pions) softer, which re- 
sults in the softening of their final transverse momentum 
spectra. At the same time, non-equilibrium corrections 
make the proton yield to increase slightly as compared to 
the case of C/s = 0.04, without changing the proton av- 
erage momentum. In the following we set Tf = 140 MeV 
and C/s = 0.08 as the optimal parameters for the freeze- 
out conditions that reproduce particle spectra and inter- 
ferometry radii in central collisions. We note, that the 
number of pions and kaons with very soft momenta is 
underestimated, while the proton yield is slightly overes- 
timated. 

The pion, kaon, and proton spectra without the bulk 
viscosity corrections at freeze-out (Eq. 4) are shown by 
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FIG. 3. (Color online) Interferometry radii in Pb-Pb collisions 
with centralities — 5% as functions of the pion pair mo- 
mentum, obtained in the viscous hydrodynamic model with 
C/s = 0.04 and Tf = 150 MeV (dashed lines), C/s = 0.04 and 
Tf = 140 MeV (dotted line), C/s = 0.08 and Tf = 140 MeV 
(solid lines). Symbols represent preliminary data of the AL- 
ICE CoUaboration [2]. 



the dashed-dotted lines in Fig. 2. Bulk viscosity correc- 
tions make the spectra softer, especially for light parti- 
cles, and introduce a correction in particle abundances, 
increasing the proton number. The corrections are im- 
portant for pions with high momenta {p± > 1.5GeV). 
They make the proton number to increase by 35%. If the 
corrections are large, formally a better ansatz for the dis- 
tribution function with bulk viscosity corrections would 
be an exponential function [37] , but the final spectra are 
very similar as when using Eq. 4. The shift in the effec- 
tive chemical equilibrium temperature is approximately 
the same for all the particles, depending only on the lo- 
cal expansion rate. A more elaborate ansatz is possible, 
with different bulk viscosity corrections for mesons, bari- 
ons, or strange particles [38], using additional parameters 
the measured particle ratios could be better reproduced. 

The spectra of identified particles in semi central colli- 
sions are well described by the hydrodynamic model (Fig. 
4). For momenta px > 1.5 GeV the pion spectra are 
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FIG. 4. (Color online) Transverse momentum spectra of vr"'", 
-ft"*", p (from to to bottom) in Pb-Pb collisions with centrali- 
ties 10 — 20% (panel a) and 30 — 40% (panel b), obtained in the 
viscous hydrodynamic model (solid lines) . Symbols represent 
preliminary data of the ALICE Collaboration [3]. 



underestimated by the hydrodynamic model with statis- 
tical particle emission. The discrepancy increases with 
centrality, and is visible for kaons as well for centrali- 
ties 30 — 40%. This effect indicates that a nonthermal 
component in the particle emission is present, e.g. jet 
fragmentation. A similar underestimation of the exper- 
imental particle yields at high momenta by the hydro- 
dynamic model is seen in peripheral Au-Au collisions at 
lower energies [37]. The proton multiplicity and spectra 
are well described by the model for different centrality 
classes. 

The production rate of strange baryons with higher 
masses S and fl is very sensitive to the chemical freeze- 
out temperature. In Fig. 5 are shown the transverse 
momentum spectra of A°, S^ and fl^ particles for two 
centrality classes. The results of the viscous hydrody- 
namic model are in qualitative agreement with the pre- 
liminary data of the ALICE Collaboration. Nonequilib- 
rium corrections in the expanding fireball increase the ef- 
fective chemical freeze-out temperature. However, the ef- 
fect is not strong enough to reproduce to observed yields 
of heavy barions. It is instructive to look at the ra- 



FIG. 5. (Color online) Transverse momentum spectra of A'' 
(dashed line), H~ (dotted line and stars), Q~ (solid line and 
diamonds) in Pb-Pb collisions with centralities — 20% (panel 
a) and 20 — 40% (panel b). The lines and the symbols rep- 
resent the results of the viscous hydrodynamic model and 
the preliminary data of the ALICE Collaboration [16] respec- 
tively. 
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FIG. 6. (Color online) Ratios of particle yields in Pb-Pb colli- 
sions with centralities 0-20%, preliminary data of the ALICE 
Collaboration (squares) [16] are compared to results of viscous 
hydrodynamic calculations (lines). 
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FIG. 7. Rapidity distributions of tt"'" (solid line), K"*" (dashed 
line) and p (dotted line) emitted in Pb-Pb collisions with 
centralities — 5%, calculated in the (3 + 1)-D viscous hydro- 
dynamic model. 



tios of pt integrated particle yields. The hydrodynamic 
model reproduces qualitatively the observed particles ra- 
tios (Fig. 6), although the nonequilibrium effects are 
described with only one parameter, the bulk viscosity 
coefficient. Deviations from chemical equilibrium at the 
freeze-out temperature Tf = 140 MeV shift the ratios of 
heavy particle yields to pion yields up. The ratio K/tt 
is very well reproduced, the ratio p/i: is overpredicted 
by 17% and the ratios S/tt, JI/tt are underpredicted by 
30 — 40%. The ratios obtained from the hydrodynamic 
simulation with Tf = 140MeV correspond aproximately 
to a chemical equilibrium temperature of 150MeV (Fig. 
6). The deviations of the calculated particles ratios from 
the observations indicate that the data cannot be de- 
scribed using a single chemical freeze-out temperature, 
as mentioned above, the agreement could be improved 
using different equilibration rates for different particle 
species. 

The (3 -I- 1)-D hydrodynamic model does not assume 
boost-invariance and gives predictions on the rapidity 
dependence of particle spectra. In practice, the initial 
density in the longitudinal direction is adjusted to repro- 
duce the final pseudorapidity distribution (Fig. 1). Once 
the transverse momentum spectra and the pseudorapid- 
ity distributions are found to be in agreement with the 
experimental observations, the rapidity distributions of 
identified particles can be reliably estimated. The rapid- 
ity distributions for pions, kaons, and protons are shown 
in Fig. 7. We find that the distributions are not boost- 
invariant. However, in the range \y\ < 1 the dependence 
of the particle densities on rapidity is very weak. This 
indicates that (2 -I- 1)-D hydrodynamic models represent 
a good approximation for the dynamics and the mech- 
anism of particle production in Pb-Pb collisions at the 
LHC at midrapidity. 

Using the hydrodynamic model with optical Glauber 
initial conditions, one can calculate the elliptic flow co- 



FIG. 8. Elliptic flow coefficient of charged particles as func- 
tion of transverse momentum for three centrality classes cal- 
culated in the viscous hydrodynamic model (lines), the sym- 
bols represent the experimental results of the ALICE Collab- 
oration [46]. 
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FIG. 9. Elliptic flow coefficient of identified particles as func- 
tion of transverse momentum calculated in the viscous hydro- 
dynamic model (lines), the symbols represent the preliminary 
results of the ALICE Collaboration [47], centrality 40 - 50%. 



efficient. The elliptic flow coefficient for charged particle 
at different centralities is presented in Fig. 8. The cal- 
culation describes well the data for soft momenta. The 
results are consistent with predictions of other viscous hy- 
drodynamic codes [6, 8-14]. Calculations of the elliptic 
flow for central collisions or of the triangular flow require 
event-by-event simulations including fluctuations in the 
initial state, and are outside the scope of this paper. An 
observable related to the freeze-out conditions for spe- 
cific particles is the elliptic flow of identified particles. 
As can be seen in Fig. 9, the hydrodynamic model with 
bulk viscosity corrections and freeze-out at 140MeV gives 
a slightly too small splitting between the elliptic flow of 
pions, kaons, and protons. 



IV. DISCUSSION 



We present (3 + 1)-D viscous hydrodynamic calcula- 
tions of the spectra of particles emitted in Pb-Pb colli- 
sions at y/sNN — 2.76 TeV. We find that the transverse 
momentum spectra are sensitive to bulk viscosity effects. 
Bulk viscosity in an exploding fireball leads to deviations 
from local equilibrium in the fluid elements. With the 
bulk viscosity coefficient (/s — 0.08 and a freeze-out 
temperature Tf — 140 MeV we find a satisfactory agree- 
ment with preliminary experimental data for the trans- 
verse momentum spectra of pions, kaons and protons. At 
the same time, the size of the fireball and the amount of 
the collective transverse flow of the fluid accumulated in 
the hydrodynamic phase are compatible with the exper- 
imental measurements of the momentum dependence of 
the interferometry radii Rout, Rside, and Riong- 

Transverse momentum spectra of identified particles in 
semi-central collisions are reproduced as well. In partic- 
ular, we find a strong transverse push visible in the px 
spectra of protons. On the other hand, bulk viscosity ef- 
fects make the pion spectra softer, in agreement with the 
experiment. Using the same freeze-out conditions we cal- 



culate the S~ and fl~ spectra and compare to the prelim- 
inary data of the ALICE Collaboration for two different 
centrality classes. Nonequilibrium corrections make the 
effective chemical freeze-out temperature to be 150 MeV 
and not 140 MeV as given by the energy density at freeze- 
out. The ratios of particle yields K/tt are very well re- 
produced, while the ratio p/ir is slightly overpredicted, 
and the strange barion abundances are underpredicted. 
This shows that an effective single chemical freeze-out 
temperature (generated dynamically in viscous hydrody- 
namics) cannot describe all the measured particle ratios. 
The elliptic flow of identified particles from the hydrody- 
namic model shows a splitting according to the particle 
mass, but slightly smaller than observed experimentally. 
Distributions of identifled particles in rapidity show that 
the system is not boost-invariant. However, in the lim- 
ited interval of central rapidities \y\ < 1 an approximate 
plateau is seen in the rapidity distributions. 
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